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ABSTRACT 

A technique is developed for the measurement of mean velocity 
direction and magnitude and the turbulence quantities in incompressible 
fluid flow within a plane, when the flow direction is known only approxi- 
mately within the plane. 

A low degree polynomial is fit to the calibration data by least 
squares, in the form: square root of velocity versus square of voltage. 
Thereby known approximate adherance to King’s law is used to advantage, 
yet accounting for small deviations from the law. 

A normal-component cooling model accounts for probe directional 
sensitivity. Calibration is performed at flow directions near those 
expected in the measurement situation, and the model is used to find a 
calibrated wire reference angle. Data taken at several flow angles 
revealed that the normal cooling model worked well in reducing the cali- 
bration data to a single curve. 

Measurements of turbulent shear stress in fully developed pipe flow 
gave excellent agreement with a theoretical prediction of total shear 
distribution. The techniques developed were also used to measure the 
velocity field of an ejector. 



MEASUREMENT OF FLOW MAGNITUDE AND DIRECTION 
BY HOT WIRE ANEMOMETER 
INTRODUCTION 

A technique is developed for the measurement of mean velocity 
direction and magnitude and the turbulence quantities in incompressible 
fluid flow within a plane, when the flow direction is known only approxi- 
mately within the plane. This situation may arise, for example, in the 
study of the development region of swirl-free axisymmetric pipe flow or 
in the study of the velocity field of a jet flow* 

The development of the method is restricted to a low level of turbu- 
lence such that the anemometer response to velocity perturbations is 
linear, i. e. , for "small” velocity perturbations. Use is made of the 
knowledge of approximate King's law behavior and of an assumed normal- 
component cooling effect. 

The technique is checked by comparing measured results with 
known behavior in fully developed pipe How. This case does constitute 
a significant simplification of procedure in that the flow direction is 
known to be parallel to the axis of the duct. Data were also taken in the 
mixing region of a ducted axisymmetric jet, where the flow direction is 
not known; and some sample results are presented. 

ANALYSIS 
Quantities Measured 

The flow quantities to be measured with respect to reference axes 
in a plane are mean velocity magnitude, U D , mean direction of flow 5, 

f~ 2 Y/2 R 

longitudinal turbulence velocity \ Ur, , transverse turbulence veloc- 

f-2 W 2 ^ Y 

ity \ v^ ? , and Reynolds stress u^ v R . 



2 


Fundamentals 

King's law for an anemometer may be written 

E 2 - E o " U eH 

where E is the voltage of the wire, and E Q is the voltage at U e ££ = 0, 
and U eff is the effective cooling velocity. This relation is developed 
from considerations of the wire energy balance including convective cool- 
ing; it is possible to evaluate the proportionality from the analysis. 

However, in practice a calibration of the wire is developed by measuring 

2 1/2 

E in a known flow. Since the relation of E versus U ' is nearly a 
linear one, the calibration data may be accurately curve -fit by a poly- 
nomial of low degree, e. g. , 

E c = 1 ( U J /2 ) - a o + h V l /2 + a 2 U c + a 3^ /2 (1) 

or 


U c /2 = s( E c) = b 0 + b l E c - b 2 E c + b 3 E c (2) 

These calibration curves may then be employed in determining wire 

sensitivities to turbulence fluctuations and in determining mean velocities. 

Use of a Normal-Sensor Probe 

Assuming that the turbulence fluctuations are "small" and that the 
calibration relation, developed in a low turbulence flow, applies instan- 
taneously in the flow, one can develop working relations. For a sensor 
which is oriented normal to the mean velocity, Fig. 1, 

U = f n (E) (3) 

and 


2 ,-,2 2 

e = sy (4) 

where S u = 9E c /3U c |y . Thus once the wire has been calibrated 
(E versus U determined) the values of U and u may be determined 

^ C 

— 2 

by measuring E and e . 
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Use of an Inclined-Sensor Probe 

For a probe lying in the xy plane and inclined to the mean velocity 
(assuming x corresponds to the direction of the mean velocity U), 

Fig. 2, the following relations apply: 


E = f (U, cp) 


e = S u + S v 
u v 


S u = 


SE 

8U 


s =iM 

Y v dcp 


u 


u 


*<p 


The relation f is fixed by calibration and may be expressed 


E c = f < U c’ <Pj 


(5) 

( 6 ) 

(7) 

(8) 

(9) 


In the general situation the probe is held so that the sensor has a fixed 
orientation with respect to reference coordinates Xp, yp, z^; these co- 
ordinates are set by the geometry of the situation (see fig 3) and may cor- 
respond, for example, to the axial, radial, and tangential directions in a 
cylindrical duct. For the present purpose the x, y and the Xp, yp planes 
are coincident, however the angle <p is unknown. In the measurement 

— 2~ p _ 

situation the unknowns are U, cp, u , v , and uv, with respect to a set of 
flow coordinates arranged such that the direction of mean velocity vector 
and the x coordinate direction coincide. It is ultimately desired to know 
the velocity and turbulence quantities relative to the coordinates Xp, yp, 
and Zp which are fixed by the geometry. Therefore 5, the angle between 
U(x) andxp, must be found. The unknowns are Up, 6, Up, Vp, and 
UpVp and at least 5 measurements are required to determine them. 

To see how these measurements are accomplished, consider the 
sensor arranged at two different orientations relative to the flow, in the 
plane of interest, Fig. 3. The angle <p Q is fixed by geometric constraints, 
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i. e. , the sensor may be held by a probe body lying along the y^ axis; 
and since cp is thereby fixed, 6 becomes the variable of interest in 
the calibration and use of the probe. 

The calibration relations may be stated: 

E cl = f i(U c ,5) (10) 

E c2 =f 2 (D c ,5) <U> 

For purposes of determining U and 6, these can be written: 

U c = g 1 (E cl ), 6 parameter (12) 

U c = g 2 ( E c2^’ 6 P arameter (13) 

These relations are established by setting up a table of E , U c , 5 in 
a flow whose direction and magnitude are controlled. However, as it is, 
the table is not very helpful in determining U and 6 when and E 2 
are measured in a flow of unknown velocity. 

The use of the calibration data would be facilitated if the angular 
dependence could be handled in equation form. It is well known that a hot 
wire is cooled primarily by the component of velocity normal to itself. 

By calibrating the probe at values of cp, i.e„, 6, near the values expected 

in the measurement situation, effects due to the velocity component along 
the wire may be accounted for in the mean. Therefore 

U eff = U sin cp 

= U sin (cp Q ± 6) (14) 

with the plus used in orientation 1 and the minus used in orientation 2. 

With the above assumption, one may write 

(U sin <p} 1/2 = g(E 2 ) = b Q + b x E 2 .+ b 2 E 4 + ( 15 ) 

and 

E 2 = f (U c sin cp} 1 / 2 (16) 



5 


The above formulations have been checked by calibrating the wire in 
orientation one in a known flow, perturbing its position by amount 5 
about angle <p Q . A sample of the result may be seen in Fig. 4. For the 
range of 6 r s chosen the data do indeed collapse to a single curve. The 
same procedure must be carried out for orientation two (no figure shown); 
the results will be the calibration relations 



These relations are then used, once and are measured in the 
measurement situation, to determine U^, V^. Using appropriate 
trigonometric identities 


2fe2\ 2(~2\ 

SilEJ + gJEj 
U R = U cos 8 - iV 17 

K 2 sin cp 

^0 

(19) 

2 cos cp Q 

(20) 

Therefore 

5 „ tan' 1 — 
U R 

(21) 

<e 1 = v 0 + p, 

(22) 

<p 2 = <p 0 - « 

(23) 

Thereby the mean velocities are determined, assuming cp Q is known. 

One may assume (p = 45°, as probe manufacturers attempt to accom 
plish, or one may determine cp Q by using the assumed normal-component 
cooling behavior: 
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[u sin {cp Q + 6)] 1/2 = g(E 2 ) (24) 

In the calibration rig U and 6 are known and E is measured. In orien 
tation one or two, one may set two known values of 5(§ A and 6g) and 


control U in the two cases so that 


e a = e b 


(25) 


For the case 6^ = 0, 6g # 0, using an appropriate trigonometric 
identity: 

U A sin cp Q = g 2 (E A ) = g 2 (E 2 ) = Ug (sin cp Q cos 5g + cos cp Q sin 6g) 
thus (^6) 


tan a = 
^o 


U B sin 6g 
U A - U B C0S 5 B 


(27) 


In practice these steps may be carried out several times, and an average 
value for cp Q may be taken. 

It remains yet to determine the turbulence quantities. Equation 
leads to the following for orientations one and two: 


e l ~ S ul u2 + 2 S ul S vl uv + S vl y2 


9 n „ „ — _2 ~~Z 

e 2 = S u2 u + 2 S u2 S v2 uv + S v2 v 


(28) 

(29) 


which may be solved for uv and v 


uv = 


2 -2 2 ^ g S 


2 2 

v2 C 1 " °vl c 2 " V°ul °v2 “ S vl S u2 


S“« e“ - SZi e 


:)u S 


VI “v2 (“ul S v2 ' S vl S u2) 


2 S.., S., 0 fs 


(30) 



( 31 ) 
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2 

These relations, however, require that u be given; this requires, 

— 9 

an additional measurement. To determine u" 1 a separate wire may be 

placed in the flow at the same location such that U is always normal to 

~~2 

the wire, thus sensitive to u . This implies aligning the wire along the 

z axis; the procedures for a normal wire may then be applied. 

Finally the turbulence components with respect to the reference 

coordinates are: 

2 2 2 2 2 

Ur = u cos 6 + v sin 6 - 2 uv cos 6 sin 6 (32) 

~2 ~2 2 ~2 2 

Vr = u sin 6 + v cos 6 + 2 uv cos 5 sin 5 (33) 

2 2 2 2 

Ur v r =(u - v ) cos 6 sin 5 + (cos 5 - sin 6) uv (34) 


The wire sensitivities for each of orientations one and two are deter 
mined as follows: 

E l = f(arg) (35) 

arg = (U c sin cp ) 1//2 (36) 

Here f represents the curve-fit of the calibration data 

2 2 3 

E c = f(arg) = a + a^ arg + a 2 arg + a^ arg + (37) 

From these 


S = 

u 


S v = 


0E 

0U 

1 0E 
U dcp 

c 


df 

U, cp Vdarg 


sin cp 


4 E (U sin cp) 
df cos cp 


1/2 


U, cp, E 


{— 

(darg 


U, cp L dar g 4 E (U sin cp) 


1/2J 


(38) 


(39) 


u, <p 9 E 


where df/darg is found by differentiating the least squares curve -fit of 


the calibration data. 
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RESULTS 

In order to evaluate the techniques which have been developed, 
measurements of uv were made in fully developed pipe flow and the 
results were compared with a theroetical prediction of the total shear 
stress distribution for such a flow. Recall that in turbulent flow the 
turbulent shear stress uv is virtually equal to the total shear stress 
everywhere in the pipe except near the wall where a laminar viscous 
shear layer exists. While the fully developed flow (5 = 0) permitted 
significant simplifications, the results provide a reasonable check of 
the techniques. 

As seen in Fig. 5 the uv results are in good agreement with the 
theoretical prediction. In the figure r represents the shear stress, 
p represents the air density, and U T represents the wall shear velocity 
found from measurement of pressure loss per unit length in the pipe, 
y is the distance from the wall and R is the pipe radius. 

Additionally Fig. 4 has already shown the effectiveness (in the 6 
range considered) of collapsing the calibration curves via the assumed 
normal cooling model. 

Results of Measurements in a Ducted Jet 

Some sample results of mean velocity magnitude and direction and 
turbulent stress in the mixing region of an axisymmetric ducted air 
jet are presented in Fig. 6 and 7. The reported variables are defined 
relative to cylindrical coordinates x, r, and 9 , which correspond to 
Xr, y^, and respectively from the section on the use of an inclined- 
sensor probe. Thus U-^ is the component of mean velocity in the axial 
direction; delta is the angle between the mean velocity vector and the 
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axial direction (positive for outward directed vector); and u R v R is the 
turbulent Reynold shear stress corresponding to axial-radial coordinates. 

The results shown are characteristic of jet mixing, with a zone of 
high radial mean velocity gradient between the inner and outer flows. 

The flow direction plot shows the outward spread of the central high 
speed flow and the inward entrainment of the outer flow into the mixing 
zone. Finally the plot of u R v R , Fig. 7, shows the high turbulent shear 
stress resulting from the large velocity gradient between the inner and 
outer flows and the very low shear stress in the potential outer flow. 

These results have been selected to demonstrate the practical 
application of the hot wire techniques. For a complete report of the 
ducted jet results, see Ref. 1. 

SUMMARY OF RESULTS 

Techniques have been presented for the measurement of mean 
velocity magnitude and direction and of turbulence quantities. The 
results were: 

1. Good agreement was found between the data and a theoretical 
prediction of shear stress in fully developed pipe flow. The pipe flow 
test did comprise a significant simplification in that the mean flow direc- 
tion was known throughout. 

2. The use of the normal-component cooling model was effective in 
accounting for the directional sensitivity of the probe in the range of 
angles considered. The model was also used to establish the reference 
angle of the probe. 



10 


CONCLUDING REMARKS 

A procedure has been developed for determining, within a plane, 
time mean velocity magnitude and direction and turbulence quantities. 

It is required that the flow direction be approximately known within this 
plane in order to effectively calibrate the probes. 

Two probes are required for the measurements: one which can be 
placed normal to the plane at all points of interest, and a second one 
which is placed in the plane at first one and then a second inclination 
relative to the flow. 

The limitations of the method are that: 

1. The probe configurations must be achievable, and obviously 
one must be able to reach the "same point” with the two probes. 

2. The flow field must be steady in the mean, since the measure- 
ments are not taken simultaneously. 

3. Due to the assumption of linear voltage response to perturba- 
tions of velocity, velocity fluctuations must be ' ’small”. 

4. In the region of interest, the sine law cooling assumption must 
adequately describe the directional effect of the flow, 

5. The calibration of the probe in a low turbulence flow applies 
instantaneously in the flow. 

6. All second order effects (e.g. , wire length) have been neglected. 

Some notable benefits of the method are that: 

1. Subject to the indicated restriction, the directional dependence 
of the calibration is handled in equation form. 

2. A digital computer may be readily applied for calibration and 
data reduction, and some on-line arrangement could be possible. 
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3. The curve fit is not forced to represent the calibration data with a 
single exponent, which would depend on the range of velocity. 

4. The reference angle of the probe (p Q is determined experimentally, 
consistent with the assumed normal component cooling model. 
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Y/R WALL DISTANCE, DIMENSIONLESS 

Figure 5. - Shear stress as a function of distance to 
wall in fully developed pipe flow. 
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RADIAL POSITION, 2y R /D ID IS DUCT DIAMETER) 

Figure 6. ~ Mean velocity magnitude and direction 
as a function of radial position in a ducted jet at 
Xd/D * 0. 5. 


RADIAL POSITION, 2R/D(D IS DUCT DIAMETER) 

Figure 7. - Shear stress as a function of radial position 
in a ducted jet at x R /D = 0. 5. 



